Transforming nonlocality into frequency dependence: a shortcut to spectroscopy 
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Measurable spectra are theoretically very often derived from complicated many-body Green's 
functions. In this way, one calculates much more information than actually needed. Here we 
present an in principle exact approach to construct effective potentials and kernels for the direct 
calculation of electronic spectra. In particular, the potential that yields the spectral function needed 
to describe photoemission turns out to be dynamical but local and real. As example we illustrate this 
"photoemission potential" for sodium and aluminium, modelled as homogeneous electron gas, and 
discuss in particular its frequency dependence stemming from the nonlocality of the corresponding 
self-energy. We also show that our approach leads to a very short derivation of a kernel that is 
known to well describe absorption and energy-loss spectra of a wide range of materials. 



PACS numbers: 71.10.-w, 71.15.Qe, 78.20.Bh 

The calculation of electronic excitations and spectra is 
one of the major challenges of today's condensed mat- 
ter physics. In fact, while Density-Functional Theory 
(DFT), especially in simple and very efhcient approxi- 
mations like the Local-Density Approximation (LDA) [l[ , 
has led to a breakthrough concerning the simulation of 
ground-state properties, the determination of electronic 
excited states is still quite cumbersome. Static ground- 
state DFT is not the adequate approach to spectroscopy. 
On the other hand, the full solution of the many-body 
Schrodinger equation for systems of more than some 
tens of atoms is today out of reach. Clearly, the de- 
scription of spectroscopy calls for the definition of suit- 
able fundamental quantities (beyond the static ground- 
state density DFT is built on), and the derivation of 
new approximations. For neutral excitations (as mea- 
sured e.g. in absorption or electron energy-loss spec- 
troscopies) one can in principle work with the time- 
dependent density, and then try to approximate this 
Time-Dependent DFT (TDDFT), for example in the adi- 
abatic LDA (TDLDA). Also Many-Body Perturbation 
Theory (MBPT) approaches [2], like the GW approxi- 
mation [31] for the one-particle Green's function or the 
Bethe-Salpeter equation (BSE) for the determination of 
neutral excitations, are today standard methods for first- 
principles calculations of electronic excitations [^. One 
of the main reasons for their success is their intuitive 
physical picture that allows one to find working strategies 
and approximations. Therefore, Green's functions ap- 
proaches are often used even when one could in principle 
resort to more simplified methods, in particular density- 
functional based ones - for example, when TDLDA breaks 
down. It is hence quite logical to try to link both frame- 
works [5] , and derive better approximations for the sim- 
pler (e.g. density-functional) approaches from the work- 



ing approximations of the more complex (e.g. MBPT) 
ones. Two prominent examples for such a procedure are 
the optimized effective potential (OEP) method ^ where 
a local and static KS potential is obtained from an in 
principle nonlocal and frequency dependent self-energy, 
and recent derivations in TDDFT of a two-point linear 
response xc kernel from the four-point BSE [7]. These 
constitute, at least at first sight, distinct methods that 
solve particular questions and leave other, very impor- 
tant points unsolved. The best example for such an open 
problem is probably the so-called "band-gap problem" [§| 
of the Kohn-Sham (KS) fH approach: one would like to 
have a simple potential that yields not only the correct 
density, but also band gaps in agreement with experi- 
ment. Although OEP in the exact-exchange approxima- 
tion (EXX) [9] had yielded some promising results in this 
direction, it is clear that a systematic solution has not yet 
been found [lOj. 

In the present work we propose a general framework 
for the definition of reduced potentials or kernels, de- 
signed for obtaining certain quantities that can other- 
wise be calculated from some A^-particle Green's func- 
tion. We show that both the abovementioned OEP and 
the TDDFT kernel from MBPT are particular cases of 
our approach. As a new level, we explore the question 
of a local and frequency dependent real potential that 
allows one to find the correct trace of the one-particle 
spectral function and therefore access direct and inverse 
photoemission (including the bandgap). We present re- 
sults of this new approach for simple metals and a model 
insulator, and discuss the relation to dynamical mean- 
field theory (DMFT) [ll| . Of particular interest is the 
conversion of nonlocality into frequency dependence that 
one encounters when an effective potential or kernel with 
a reduced number of spatial degrees of freedom is used 
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to describe excitations. 

In this general scheme we now suppose that we want 
to calculate the quantity T that is a part of the informa- 
tion carried by the iV-particle Green's function G. We 
symbolically express this relation as: T = p{G}. In the 
following we will specify the formula for the one-particle 
G, whereas the completely analogous case of two-particle 
Green's functions will be briefly discussed at the end. 
As a well-known example we take the electronic density 
T = p for which the "part" to be taken is the diagonal of 
the one-particle G: p{r) = —iG{r,r,t,t^). We then in- 
troduce another Green's function Gt which has the part 
p{.} in common with G. T = p{Gt}. We also suppose 
that Gt is associated to an effective potential Vt accord- 
ing to Gt = (j^ — Ho — Vt)~'^, where Hq is the Hartree 
Hamiltonian [2l|. The full Green's fimction G and the 
new Gt are linked by a Dyson equation: 



G = Gt + Gt{^- Vt) G. 



(1) 



We now take the part of interest p{.} of this Dyson equa- 
tion. This yields the condition: 



p{Gt (S - Vt) G} = 0. 



(2) 



The aim is now to make an ansatz for Vt with a sim- 
pler structure than S, and for which ([2]) can be solved. 
For the example where T is the density p, a static and 
local potential can do the job. In fact in this example 
Vt is the exchange-correlation potential Vxcir), Gt is 
the Kohn-Sham Green function Gks, and ([2]) is a well- 
known result, derived by Sham and Schliiter in [12]. It 
has subsequently been extended to time-dependent ex- 
ternal potentials [l3| . and einployed in many different 
contexts (see e.g. Refs. [i][l3]). The Sham-Schliiter 
equation (SSE) is still implicit since self-consistent in Vxc 
and Gks- Most often it is hnearized setting G = Gks 
everywhere, including in the construction of S. This so 
called linearized Sham-Schliiter equation, that can also 
be derived using a variational princij Die [Hill, is the 
central equation of the OEP method [6|. In particular, if 
one uses for S the exchange-only approximation S = S^; , 
where Ti^ — iGksv {v is the bare Coulomb interaction), 
one obtains the so called exact-exchange approximation 
to Vxc • Obviously, by construction one obtains in this 
way a good description of the density, whereas there is 
no reason for other features of the Green's function, in 
particular the band gap, to be correct. However, in our 
much more general scheme one can easily go beyond the 
case of solely the density p{r). According to the choice 
o{p{.}, the ansatz ioT Vt has then to be modified. Below 



we will first explore the problem of electron addition and 
removal. 

Electronic structure, defined as electron addition and 
removal energies, is measured by experiments like direct 
or inverse photoemission. To first approximation Q these 
experiments measure the trace of the spectral function 
A{ri,r2,uj) = ^\lmG{ri,r2,uj)\. In other words, for the 
interpretation of photoemission spectra one doesn't need 
the knowledge of the whole Green's function G, but just 
the imaginary part of its trace over spatial coordinates to- 
gether with its full frequency dependence. We will hence 
look for a potential that is simpler than the full self- 
energy but yields the correct trace of the spectral func- 
tion. It is reasonable to add the condition that also the 
density should be correct, which means that the diagonal 
in real space, and not only its integral, is fixed. Following 
the general scheme, we introduce a new Green's function 
GsF ~ — Ho — Vsf)~^ stemming from a potential Vsf 
such that: lmGsF{ri,ri,uj) = ImG(ri, ri, w). What de- 
grees of freedom are needed in V5F? A natural assump- 
tion is that J^^F should be local in space, but frequency 
dependent [l7|[18|. It is also possible to choose Vsf to 
be real. With this ansatz Eq. ([2|) yields: 

VsF{ri,uj) = J dr2dr3dri^^^{ri,r4,uj) x 

lui{GsFir4,r2,uj)T.{r2,r3,u)G{r3,r4,u)} (3) 

where ^(ri,r2,w) = lm{GsF{ri,r2,uj)G{r2,ri,uj)}. The 
solution exists if ^ is invertible. Whether this is always 
the case is a delicate question and beyond the scope of the 
present work. Equation ([3]) shows that Vs f should indeed 
be frequency dependent unless S is static and local (in 
that case the w-dependent terms cancel trivially). Thus, 
in general a local static (KS) potential will not be able to 
reproduce the spectral function, whereas T/gi?(r, w) is the 
local potential that will yield the correct bandgap and 
the correct density p{r) of the system. 

At this point it is interesting to compare our con- 
struction with the approach of the Spectral Density- 
Functional Theory (SDFT) [l^, where the key vari- 
able is the short-range part of the Green's function: 
Giocir, r', oj) = G{r, r' , oj) Q{il.c), where <d{ilc) is 1 when 
r is in the unit cell and r' inside the volume flc, and is 
otherwise (see Fig. 1 of Ref. [l^). A new Green's func- 
tion GsDFT = {^^ — Hq — Vsdft)~^ can be introduced 
such that GsDFT = Gioc where G/oc is different from 
0. Using this property of Gsdft, we find that the (in 
general complex) potential VsdfTi defined in the volume 
nc, is: 



V, 



SDFT 



ir5,re,u) 



dridr2 / dr^dr^G 



SDFT 



{r5,ri,uj)GsDFT{ri,r3,uj)j:{r3,r4,uj)G{r4,r2,Uj)G ^{r2,rQ,uj) (4) 
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where G~^, if it exists, is the local inverse of G in flc 
(while would be the full inverse, defined in the whole 
space). In principle SDFT is a formally exact theory. 
The most common approximation to SDFT is the dy- 
namical mean- field theory [T3|. In this perspective |T§| 
DMFT can be viewed as consisting in the assumption 
of taking for the interaction energy functional of SDFT 
^SDft[Gioc] the form of functional $mb[G] for the full 
Green's function G. This means that the shorter the 
range of the self-energy is, the better this approximation 
becomes. In particular, in the limit case that S is com- 
pletely localized in flc, then Vsdft and Gsdft coincide 
respectively with S and G. The interesting situation is of 
course when this is not true. In fact, our Vsf of Eq. ([3]) 
corresponds to the case where ftc ^ so that this condi- 
tion is certainly not fulfilled. Then, as we will illustrate 
below, the nonlocality of E will strongly influence Vsf 
and in particular lead to a frequency dependence which 
is not the frequency dependence of S itself. This will to a 
certain extent also be true for any Qq of finite range, so 
that the following discussions may also give useful insight 
for research in the field of DMFT. 

To illustrate the frequency dependence of Vs f we con- 
sider the case of homogeneous systems, where all local 
quantities (like Vxc and Vsf) are constant in space. In 
particular, the xc potential is Vxc = S(p = pf,^ = 0), 
while, since Vsf = Vsi?(w), one can directly write (here 
and throughout the paper we adopt atomic units): 

\lmGsFir,r,Lo)\^2 J ^nS{Lu + ^l ~ - Vsf{l^)) 

= — eiio + fi- VsFiio))V^ + M - Vsf{uj). (5) 

Requiring that this is equal to |ImG(r, r, a;)|, one finds a 
unique local potential Vsf (for Vsf{^^) < uj + fJ-): 



Vsf{^) = oj + 



/ TT 

[—^ |ImG(r, r,w) 



(6) 



It is first of all interesting to consider the case of simple 
metals, that can be modelled by a homogeneous elec- 
tron gas. Here, we will assume a static but nonlocal 
self-energy T,x{r-r') = iG{r - r' ,t - t+)vxir - r'), 
that is a screened-exchange like form where the screened 
Coulomb potential is v\{r — r') = v{r — r')e~^'^~^ For 
a larger screening length A, S is more effectively nonlo- 
cal, since less screened. A tunes the effective range of the 
interaction: from A = (E = 0, Hartree approximation), 
to A = oo (unscreened Hartree- Fock). Using the relation 
^ we calculate Vsf{'-^) for different values of A, i. e. for 
different nonlocality ranges of E (see Fig. [T]). We find 
that the spatial nonlocality of the static self-energies is 
completely transformed into the frequency dependence of 
Vsf- This is an essential property of Vsf, which radi- 
cally distinguishes Vsf from the static Vxc- In particular, 
looking at Fig. [U one observes that the more nonlocal Y,\ 
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FIG. 1: (Color online) Transformation of nonlocal statically 
screened exchange self-energy Sa(p) (upper panel) to the fre- 
quency dependent local Vsf(i^) (bottom panel), for differ- 
ent screening lengths A in HEG {\tf ~ (4pf/7I")~^'^^ is the 
Thomas- Fermi length), and densities corresponding to alu- 
minium (lines) and sodium (open symbols). The key is com- 
mon to both panels. The (^-dependence of Vsf is stronger 
for more nonlocal self-energies. Below band-edges, where 
ImG = 0, Vspii-o) has been defined continuous and equal to a 
constant > oj + ^i. In the inset: comparison between the exact 
and approximated solutions © and ^ in the Hartree-Fock 
case. 



is, the more dynamical Vsf becomes. This is confirmed 
by comparing two materials: aluminium and sodium. 

The HEG is a prototype of metallic systems. Let us 
consider the simplest model insulator introduced by Call- 
away [20[. The Callaway's model is obtained from HEG 
by inserting a gap A between the occupied "valence" 
states 01, (r) (p < pp) and the "conduction" states 0c('') 
(p > pf)- Within the Green's functions formalism this 
corresponds to a nonlocal static "scissor" self-energy 

S(r,r') = A^0,(r)</>:(/), (7) 

C 

which rigidly shifts all conduction states, and produces 
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a gap between the valence, = ep = fJ-, and the con- 
duction, Ec = sp + band edges. Due to the homo- 
geneity of the model, Eq. ^ is still applicable. From 
Eq. one gets that Vsf{^) — for lu < E^ — ^ ^ 
and 1/sf(w) = A ior uj > Ec ~ h — A. Between a; = 
and w ^ A, where ImG(cj) = 0, there are many different 
choices for Vspi^^)- from (O it is enough that Vsf{^^) > 
Lu + ^to get IrnGspi^) — 0- In any case a static constant 
potential, like Kcc cannot produce the correct bandgap, 
because at Ey — fi and Ec ~ fJ- it should assume different 
values. The only effect of Vxc with respect to Hq would 
be just a rigid shift of the whole bandstructure and no 
gap could be opened. In fact, A in ([7|) corresponds to the 
derivative discontinuity of Vxc when an electron is added 
to the system: A = V^^^^\r) - V^^\r) ^. Our sim- 
ple example demonstrates that this discontinuity, which 
constitutes the difference between the true quasiparticle 
gap and the Kohn-Sham one, is accounted for by the po- 
tential VsF rightly through its frequency dependence. 

Knowing the exact solutions, one can also verify how 
common approximations to SSE perform with these 
strongly frequency dependent functions, in particular 
when we linearize the SSE by setting G = Gsf every- 
where (in DFT this would correspond to the OEP ap- 
proach) . The linearized SSE is still self-consistent in Vsf 
and Gsf- Yet we can approximate it at the first order 
setting G and Gsf equal to the Hartree Gh- In the 
Hartree-Fock case in HEG, we get, for a; > —Ef- 
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(8) 



which shows a logarithmic divergence at the Fermi energy 
w = and hence is very different from the exact solution 
(see the inset of Fig. [T]) . The approximated Vs f in HEG 
is distant from the exact one also when the linearized SSE 
is solved self-consistently. The agreement should improve 
when screening is taken into account, since then Vsf is 
less frequency dependent and Gsf is much closer to G 
and Gu . Already for sodium the agreement is better than 
for aluminium. In any case, these results demonstrate 
that the linearization of the generalized SSE is a delicate 
procedure. In the following we will however show that 
it works very well in certain cases, for example for the 
kernel j^c of TDDFT. 

Up to now we have considered Dyson equations only 
for the one-particle Green's function. The generalization 
of SSE can however be applied to any Dyson-like equa- 
tion, for example involving two-particle Green's func- 
tions, like the four-point reducible polarizability ^X-A/b- 
In this way we can derive an exact equation for the ker- 
nel Uc = fxc + fxc of TDDFT. In particular, follow- 
ing Ref. J5i|j we concentrate on the term jxc that de- 
scribes electron-hole interactions. This means that in 
the TDDFT polarizability "^xtd the "gap opening" con- 
tribution due to the term /ic'' is already included by us- 



ing = —iGG instead of the KS independent parti- 
clc polarizability. Jxc is the two-point kernel that gives 
the correct two-point polarizability [l^l- '*Xtd(1122) = 
"'xmb(1122) = x(12). Applying this condition to the 
Dyson equation "^xmb = ^Xtd[Fmb - /i?^] "^Xmb yields 
the generalized SSE, and hence the exact expression of 

f(2). 
Jxc ■ 



/i?'(34) = x"'(31) Vd(1156)Fms(5678)x 

X V/i3(7822)x-i(24), (9) 

where Fmb = i6T,/SG is the exchange-correlation part of 
the kernel of the Bethe-Salpeter equation, which, in the 
framework of ab initio calculations [4| , is most frequently 
approximated by: FA/i3(1234) = -W'* {12)6 {13)6 {2 A) 
where W'*(12) = W{ri,r2,Lu = 0)6{t2 - t+) is the 
statically screened Coulomb interaction. An approxi- 

(2) . 

mated expression for fxc , obtained from the linearized 
generalized SSE, where we set x(12) = X?^(12) and 
xHl; 23) = 4xmb(1123) - ^Xtd{1123) - ^X^lUSS) - 
--iG(12)G(31), is then: 



-\Lo)xH^)W{0)\{u;){xrr'i^) 



(10) 

In this equation only the frequency dependence of the 
various terms has been put into evidence, since it is 
mostly interesting to note that, although W and hence 
Fmb are static, fxc is frequency dependent unless W is 
short-ranged in which case the frequency dependence of 
the other components cancels. This recalls the analogous 
transformation of static self-energies into frequency de- 
pendent potentials discussed above. Equation pO|) rep- 
resents a new derivation of a two-point linear response 
exchange-correlation kernel that has previously been ob- 
tained in several other ways Q. It has been shown to 
yield absorption or energy-loss spectra of a wide raiige 
of materials in very good agreement with experiment [7| . 
This means also that the linearization of the generalized 
SSE in this case turns out to be a very good approxi- 
mation. The present derivation is particularly quick and 
straightforward, showing one of the advantages of the 
generalized Sham-Schliiter equation formulation. 

In conclusion, in this paper we propose a shortcut for 
the calculation of electronic spectra, based on a general- 
ization of the Sham-Schliiter equation [l^ . In particular, 
we have introduced a local and real potential for pho- 
toemission, with a frequency dependence stemming both 
from the frequency dependence and from the non-locality 
of the underlying self-energy. We have illustrated some 
features at the example of sodium and aluminium. We 
have also applied the approach to the derivation of an 
exchange-correlation kernel for the calculation of absorp- 
tion spectra. This work opens the way to explore new 
approximations for simplified potentials and kernels that 
can be employed to calculate a wide range of electronic 
spectra. 
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